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We review recent developments in non-perturbative field theory using modal field 
methods. We discuss Monte Carlo results as well as a new diagonalization tech- 
nique known as the quasi-sparse eigenvector method. 

1 Introduction 

Of the many approaches to non-perturbative quantum field theory, we can 
identify two general computational strategies. The first is the method of 
Monte Carlo. The main advantages of this approach is that it can treat 
many higher dimensional field theories, requires relatively little storage, and 
can be performed with massively parallel computers. The other strategy is the 
method of explicit diagonalization. The strong points of this method are that 
it is immune to the fermion sign problem, can handle complex-valued actions, 
and yields direct information about the spectrum and eigenstate wavefunc- 
tions. 

Modal field theory is a simple Hamiltonian framework which can acco- 
modate either computational strategy. The first step is to approximate field 
theory as a finite-dimensional quantum mechanical system. The approxima- 
tion is generated by decomposing field configurations into free wave modes 
and has been explored using both spherical partial waves 1 and periodic box 
modes 2 . From there we can analyze the properties of the reduced system 
using Monte Carlo, diagonalization, or some other computational method. In 
this short review we present two different approaches using the modal field 
formalism which address the example of 4 theory in 1 + 1 dimensions. We 
begin with a summary of the Monte Carlo calculation of the critical behavior 
of <f) A theory in Euclidean space. We then discuss the computational chal- 
lenges of a diagonalization-based approach and conclude with a calculation 
of the lowest energy states of the theory using a diagonalization technique 
known as the quasi-sparse eigenvector method. 

2 Monte Carlo 

In this section we review the basic features of modal field theory in a pe- 
riodic box and discuss spontaneous symmetry breaking in Euclidean two- 
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dimensional (p 4 theory using the method of diffusion Monte Carlo. The 
field configuration <f) is subject to periodic boundary conditions (f>(t, x — L) = 
(f>{t, x + L). Expanding in terms of periodic box modes, we have 

<Kt,x) = Jj t J2 Mt)e inwx/L - (1) 

n=0,±l,... 

The sum over momentum modes is regulated by choosing some large positive 
number iV max and throwing out all high-momentum modes <f) n such that \n\ > 
AT max . In this theory renormalization can be implemented by normal ordering 
the </> 4 interaction term. After a straightforward calculation (details shown in 
2 ), we find that the counterterm Hamiltonian has the form 

E ( 2 ) 



4!2L 

n=-AT max ....Ar, 



where 



h = E *b -n = \/4^+M 2 - (3) 

We represent the canonical conjugate pair <fr n and d ^Z n using the Schrodinger 
operators q n and Then the functional integral for <f> 4 theory is equiv- 

alent to that for a quantum mechanical system with Hamiltonian 

H= E [-heh^ + ^l-^-nln] (4) 

rt=— JV max ,...JV max 



n 1 +n 2 +n 3 +n 4 =0 

The existence of a second-order phase transition in two-dimensional (j) 4 
theory has been discussed in the literature 3 . The phase transition is due 
to tfi developing a non-zero expectation value and the spontaneous breaking 
of (j) — ► — <P reflection symmetry. It is believed that this theory belongs to 
the same universality class as the two-dimensional Ising model and therefore 
shares the same critical exponents. From the Ising model we expect to find 

u = l, 0=\. (5) 

v is the exponent associated with the inverse correlation length or, equiv- 
alently, the mass of the one-particle state, m. It is defined by the behavior 
of m near the critical coupling A c , 

m~(\ c -\y. (6) 



mftqse: submitted to World Scientific on February 1, 2008 2 



We determine the behavior of the mass as we approach the critical point 
from the symmetric phase of the theory. All computations are done using 
the method of diffusion Monte Carlo (DMC). The idea of DMC is to model 
the dynamics of the imaginary-time Schrodinger equation using the diffusion 
and decay/production of simulated particles. The kinetic energy term in the 
Hamiltonian determines the diffusion rate of the simulated particles and the 
potential energy term determines the local decay/production rate. A self- 
contained introduction to DMC can be found in 4 . 

(3 is the critical exponent describing the behavior of the vacuum expecta- 
tion value. In the symmetric phase the vacuum state is unique and invariant 
under the reflection transformation 0^—0 (or cquivalently q n — ► — q n , for 
each n). In the broken-symmetry phase the vacuum is degenerate as L — > oo, 
and qo, the zero-momentum mode, develops a vacuum expectation value. In 
the L — > oo limit tunnelling between vacuum states is forbidden. One ground 
state, |0 + ) , is non-zero only for values qo > and the other, |CT), is non-zero 
only for qo < 0. We will choose |0 + ) and |0~) to be unit normalized. We 
can determine (3 from the behavior of the vacuum expectation value as we 
approach the critical coupling, 

(O+|0|O+)~(A-A C )' 3 . (7) 

In our calculations we used (L, N max ) = (2.5n, 8), (2.5w, 10), (5n, 16), and 
(5-7T, 20). For convenience we measure all quantities in units where fj, = 1. For 
each set of parameters L and iV max , the curves for m and (0 + | cj> |0 + ) near the 
critical coupling were fitted using the parameterized forms 

™ = «(¥-^f (8) 

and 

<0 + k|0 + )=&(A-^) /3 . (9) 

Combining the results for the various values of L and iV max and taking into 
account finite-size effects, we can extrapolate to the L, iV max — > oo limit. The 
results we find are 

= 2.5 ± 0.2 ± 0.1, v = 1.3 ± 0.2 ± 0.1, a = 0.43 ± 0.05 ± 0.02, (10) 
4^ = 2.5 ±0.1 ±0.1, (3 = 0.13 ±0.02 ±0.01, b = 0.71 ± 0.04 ± 0.03. 

The first error bounds include inaccuracies due to Monte Carlo statistics, 
higher energy states (for the mass curves), and extrapolation. The second 
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error bounds represent estimates of the systematic errors due to our choice 
of initial state, time step parameter, and bin sizes in the DMC simulations. 
Our results for the critical exponents are consistent with the Ising model 

predictions (5). The results for the critical coupling and ^j^- are in 
agreement with a recent lattice computation 5 

If = 2.56±;C? (11) 
as well as the discrete light-cone quantization result 6 

*f » (3+V5)* „ 2.48. (12) 

3 Diagonalization 

In principle diagonalization techniques can provide significant detailed in- 
formation about a quantum system. The basic idea is to diagonalize the 
Hamiltonian matrix defined over a truncated finite-dimensional Fock space. 
In practise however few field theories are computationally feasible and typi- 
cally only theories in 1 + 1 dimensions. Sparse matrix techniques such as the 
Lanczos or Arnoldi schemes allow us to push the dimension of Fock space to 
about 10 5 or 10 6 states. This may be sufficient to do accurate calculations for 
4>\ + i theory near the critical point ^ w 2.5 for larger values of L and iV max . 
It is, however, near the upper limit of what is possible using current computer 
technology and existing algorithms. Unfortunately field theories in 2 + 1 and 
3 + 1 dimensions will require much larger Fock spaces, probably at least 10 12 
and 10 18 states respectively. In order to tackle these larger Fock spaces it is 
necessary to venture beyond standard diagonalization approaches. 

In view of this difficulty one could ask how the Monte Carlo method avoids 
this problem. The answer is by means of importance sampling. While the 
space of all configurations is impossibly large, it is often adequate to sample 
only the most important configurations. Unfortunately Monte Carlo tech- 
niques are not effective for a number of important problems. These include 
systems with dynamical fermions and computations involving complex-valued 
actions. In general difficulties arise when the functional integral measure is 
not positive definite and produce contributions with oscillating signs or phases. 

In 7 a new importance sampling method was proposed for use within 
a diagonalization scheme. The aim was to exploit the extreme sparsity of 
the Hamiltonian matrix, a result of the restricted form that characterizes local 
renormalizable interactions. If N is the number of non-zero entries per row or 
column in the Hamiltonian matrix, it is shown in 7 that a typical eigenvector 
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is dominated by its largest y/N elements. There are some exceptions to this 
rule related to degeneracies which are discussed in 7 . The algorithm for 
diagonalizing H is as follows. We start by choosing a set of orthonormal basis 
vectors such that the Hamiltonian matrix is sparse and the eigenvectors 
are quasi-sparse. The term quasi-sparse indicates that the norm of a vector is 
dominated by the contribution from a small fraction of its components. The 
remaining steps are as follows: 

1. Select a subset of basis vectors {e^ , • • • , ei n } and call the corresponding 
subspace S. 

2. Diagonalizc H restricted to 5* and find one eigenvector v. 

3. Sort the basis components (e^ \ v) according to absolute size and throw 
away the least important basis vectors. 

4. Replace the discarded basis vectors by new basis vectors. These are 
selected at random from a pool of candidate basis vectors which are 
connected to the old basis vectors through non- vanishing matrix elements 
of H. 

5. Redefine S as the subspace spanned by the updated set of basis vectors 
and repeat steps 2 through 5. 

If the subset of basis vectors is sufficiently large, the exact eigenvectors will be 
stable fixed points of the update process. We will refer to this diagonalization 
technique as the quasi-sparse eigenvector (QSE) method. 

We now test the QSE method on (j) 4 theory in 1 + 1 dimensions. The 
theory can be treated rather well using conventional sparse matrix techniques 
and a Fock space of about 10 6 states. However the test we propose will be 
more stringent. The size of Fock space and the number of non-zero transitions 
per state for theory in 3 + 1 dimensions are the cubes of the corresponding 
numbers in 1 + 1 dimensions. Therefore a more useful test is to solve the 1 + 1 
dimensional system in a manner that can be directly generalized and expanded 
to the 3+1 dimensional problem. Setting an upper limit of 10 6 states for the 
3 + 1 dimensional case, we will attempt to solve the 1 + 1 dimensional system 
using no more than 100 (cube root of 10 6 ) basis vectors at a time. 

We again use the Hamiltonian in (4) and set L — 5n and N max = 20. 
This corresponds with a momentum cutoff scale of A = 4. We also place 
auxiliary constraints on the states in our momentum Fock space. We keep 
only those states with < 13 particles and kinetic energy < 2A. A precise 
definition of the kinetic energy cutoff is provided in 8 , and in terms of the 
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parameter K max introduced there we use K m!iX = 41. In our calculations we 
also restrict our attention to the zero momentum sector. 

With these constraints our Fock space contains about 2 x 10 6 states and 
the Hamiltonian matrix has about 10 3 transitions per state. However, we 
will restrict our QSE calculation to include only 100 basis vectors at a time. 
Results for the energy eigenvalues are shown in Fig. 1. From our Monte 
Carlo calculation we know that the theory has a phase transition at w 2.5 
corresponding with spontaneous breaking of the 0-^—0 reflection symmetry. 
In the broken phase there are two degenerate ground states and we refer to 
these as the even and odd vacuum states. In Fig. 1 we see signs of a second 
order phase transition near ^ w 2.5. Since we are working in a finite volume 
the spectrum is discrete, and we can track the energy eigenvalues as functions 
of the coupling. Crossing the phase boundary, we see that the vacuum in 
the symmetric phase becomes the even vacuum in the broken phase while the 
one-particle state in the symmetric phase becomes the odd vacuum. The 
energy difference between the states is also in agreement with the Monte 
Carlo calculation of the same quantities. The state marking the two-particle 
threshold in the symmetric phase becomes the one-particle state above the 
odd vacuum, while the state at the three-particle threshold becomes the one- 
particle state above the even vacuum. These one-particle states should be 
degenerate in the infinite volume limit. One rather unusual feature is the 
behavior of the first two-particle state above threshold in the symmetric phase. 
In the symmetric phase this state lies close to the two-particle threshold. But 
as we cross the phase boundary the state which was the two-particle threshold 
is changed into a one-particle state. Thus our two-particle state is pushed up 
even further to become a two-particle state above the even vacuum and we 
see a pronounced level crossing. 

A conventional sparse matrix algorithm would require about 10 GB RAM 
and 10 5 Gflops to determine the energy and wavefunctions of the lowest lying 
states for this system. The QSE method was able to perform the same task 
with about 30 KB RAM and 10 2 Gflops. The reductions in floating point op- 
erations and memory are about three and six orders of magnitude respectively 
and show the advantages of QSE diagonalization over other numerical tech- 
niques. Other large non-perturbative systems should now also be amenable 
to numerical solution using this method. 
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Figure 1. Energy eigenvalues of <fi\ , l as functions of the coupling constant. 
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